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Ion-acoustic modes propagating in unmagnetized dusty plasmas are studied by applying a generic 
collisionless fluid model. An Eulerian-to-Lagrangean variable transformation leads to a new system 
of evolution equations, which may be combined into a single (rather complex) equation for the mean 
ion density. The linear and weakly nonlinear oscillations regimes are studied and their modulational 
stability is investigated. This study extends known previous results on the Lagrangian description of 
electron- or ion-acoustic electrostatic plasma modes. The relevance with the description of dust-ion 
acoustic waves propagating in a dusty plasma is also discussed. 
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I. INTRODUCTION 

Dusty plasmas (DP) have recently attracted a great deal of attention due to a variety of new phenomena observed in 
them and the novel physical mechanisms involved in their description 0,0- In addition to known plasma electrostatic 
waves |2j, e.g. ion-acoustic waves (IAW), new oscillatory modes arise in DP 0,12, among which the dust-ion acoustic 
wave (DIAW), essentially corresponding to ion-acoustic (IA) oscillations in a strong presence of dust. 

The linear properties of the IAWs and DIAWs have been quite extensively studied and now appear well understood 
0. As far as nonlinear effects are concerned, various studies have pointed out the possibility of the formation of 
DIAW-related localized structures, due to a mutual compensation between nonlinearity and dispersion, including 
small-amplitude pulse solitons, shocks and vortices Furthermore, the propagation of nonlinearly modulated DIA 
wave packets was studied in Ref. [5| , in addition to the formation of localized envelope soliton-modulated waves due 
to the modulational instability of the carrier waves. A very interesting known approach is the Lagrangean description 
of a nonlinear wave profile. In the context of electrostatic plasma waves, this formalism has been employed in studies 
of electron plasma waves 0, 0, and, more recently, ion-acoustic and dust-acoustic 01 waves. Our aim here is 
to extend previous results by applying the Lagrangean formalism to the description of nonlinear DIAWs propagating 
in dusty plasmas. 

II. THE MODEL 

We shall consider the nonlinear propagation of dust-ion-acoustic waves in a collisionless plasma consisting of three 
distinct particle species 'a': an inertial species of ions (denoted by V; mass m^, charge qi = +Zie; e denotes the 
absolute of the electron charge), surrounded by an environment of thermalized electrons (mass m e , charge — e), and 
massive dust grains (mass M, charge qa = sZ^e, both assumed constant for simplicity); Zj, denotes the charge state 
of dust grains; we leave the choice of dust grain charge sign s = qd/\qd\ (= — 1/ + 1 for negative/positive dust charge) 
open in the algebra. Charge neutrality is assumed at equilibrium. 

The hydrodynamic-Poisson system of equations which describe the evolution of the ion 'fluid' in the plasma are 
considered. The ion number density rii is governed by the continuity equation 

^ + V(n lUl ) = 0, (1) 



* Proceedings of the International Conference on Plasma Physics - ICPP 2004, Nice (France), 25 - 29 Oct. 2004; contribution P3-045; 

available online at: http://hal.ccsd.cnrs.fr/ccsd-00001891/en/ . 
t On leave from: U.L.B. - Universite Libre de Bruxelles, Physique Statistique et Plasmas C. P. 231, Boulevard du Triomphe, B-1050 

Brussels, Belgium; also: Faculte des Sciences Apliquees - CP. 165/81 Physique Generale, Avenue F. D. Roosevelt 49, B-1050 Brussels, 

Belgium; 

Electronic address: ioannisatp4.rub.de 
t Electronic address: ps0tp4.rub.de 



2 



where the mean velocity obeys 



duj Zje Z^e 

— - + u, • Vu, = — E = V $ . 2 

at rrii rrij 



The electric field E = —V $ is related to the gradient of the wave potential <I>, which is obtained from Poisson's 
equation V • E = An q s n si viz. 

V 2 $ = 4tt e (n c + n h - Zi th) ■ (3) 

Alternatively, one may consider 

dE 



dt 



= -47T 2J <7a "a "a • (4) 



We assume a near-Boltzmann distribution for the electrons, i.e. n e w n Ci o exp(e$/fcsr e ) (T e is the electron temper- 
ature and fcs is Boltzmann's constant). The dust distribution is assumed stationary, i.e. rid ~ const.. The overall 
quasi-neutrality condition at equilibrium then reads 

Zi n lj0 + sZ d n d - n ej0 = . (5) 



A. Reduced Eulerian equations 



By choosing appropriate physical scales, Eqs. (l)-(3) can be cast into a reduced (dimensionless) form. Let us define 
the ion-acoustic speed c s = (ksTe/mi) 1 / 2 . An appropriate choice for the space and timescales, L and T = L/c s , are the 
effective Debye length = (^T e / '4tt Zfrii^e 2 ) 1 ^ 2 = c s /uj Pi i and the ion plasma period u>~\ — (ATrnifiZfe 2 /mi)^ 1 ' 2 , 
respectively. Alternatively, one might leave the choice of L (and thus T = L/c s ) arbitrary - following an idea suggested 
in Refs. [flLllfl^ - which leads to the appearance of a dimensionless dispersion parameter S = 1/{lo p ^T) = \d/L in the 
formulae. The specific choice of scale made above corresponds to 5 = 1 (implied everywhere in the following, unless 
otherwise stated); however, we may keep the parameter 5 to 'label' the dispersion term in the forthcoming formulae. 

For one-dimensional wave propagation along the x axis, Eqs. JTJ - J2J can now be written as 



0. 



dn d(nu) 

dt dx 
du du _ , 

where all quantities are dimensionless: n — ni/riifi, u = ui/vq and <p = ^/^o; the scaling quantities are, respectively: 
the equilibrium ion density n^o, the effective sound speed vq = c s (defined above) and $o = kBT e /(Zie). The 
(reduced) electron and dust background density h is defined as 

ft=Jlf-e^ + .^=^ + l-/ 1 , (7) 
ZiTlifi Ziiiifl 

where we have defined the DP parameter /x = n e> o/(-^i^j,o)> an d made use of Eq. (J2J). Note that both n and h reduce 
to unity at equilibrium. 

We shall define, for later reference, the function f((fi) = n - given by Eq. J7J) - and its inverse function 

f-\ x )= Zi \nU + ^pj=9{x), (8) 

viz. f(4>) — x implies (f> = f^ 1 (x) = g{x). 

We note that the dependence on the charge sign s is now incorporated in fj, = 1 + sZ^n^Q) '(Zjn^o); retain that 
fx < 1 (/x > 1) corresponds to negative (positive) dust. Remarkably, since the dust-free limit is recovered for fi = 1, the 
results to be obtained in the following are also straightforward valid for ion-acoustic waves propagating in (dust-free) 
e-i plasma, upon setting ji = 1 in the formulae. 
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The well-known DIAW dispersion relation u 2 — c 2 k 2 / {k 2 \ 2 D + 1) is obtained from Eqs. to J5J. On the other 
hand, the system yields the reduced relation co 2 = k 2 /(S 2 k 2 + 1), which of course immediately recovers the former 
dispersion relation upon restoring dimensions (regardless, in fact, of one's choice of space scale L; cf. definition of 
S). However, some extra qualitative information is admittedly hidden in the latter (dimensionless) relation. Should 
one consider a very long space scale L » An (i.e. S -C 1), one readily obtains lu ~ k (unveiling the role of S as a 
characteristic dispersion control parameter). Finally, the opposite limit of short L (or infinite 6) corresponds to ion 
plasma oscillations (viz. ui = lo p ^ = constant). 

B. Lagrangean description 

Let us introduce the Lagrangean variables {£, r}, which are related to the Eulerian ones {x, t} via 

£ = x - f u(^T')dr' , r = t. (9) 



Jo 

See that they coincide at t — 0. Accordingly, the space and time gradients are transformed as 

d/dx -^a^d/dS,, d/dt ^d/dr-a^u d/d£ , 
where we have defined the quantity 

Br, [ T B 

a^,r) = — =l + ^ dr'-u(tr'). (10) 

Note that the convective derivative D = d/dt + u d/Bx is now plainly identified to d/dr. Also notice that a satisfies 
a(£, r = 0) = and 

da(£,T) _ du(£,r) 

dr d£ [ ' 

As a matter of fact, the Lagrangean transformation defined here reduces to a Galilean transformation if one suppresses 
the evolution of u, i.e. for u — const, (or du/dr = du/d^ — 0, hence a = 1). Furthermore, if one also suppresses the 
dependence in time r, this transformation is reminiscent of the travelling wave ansatz /(x, t) — f(x — vt = s), which 
is widely used in the Sagdeev potential formalism 0, [n| . 

The Lagrangean variable transformation defined above leads to a new set of reduced equations 

»(£,r) = ar 1 ^, t)<,0) (12) 

= -^a-^,r)^fl (13) 
mi d£ 

<*~Ht>T) dE[ ^ T) = 47rZ i e[n(C ! r)-nn i , ] (14) 

— -a^u—jE^T) = -4*cZien(Z,T)u(£,T), (15) 

where we have temporarily restored dimensions for physical transparency; recall that the (dimensionless) quantity h, 
which is in fact a function of <f>, is given by |7f . One immediately recognizes the role of the (inverse of the) function 
cx(£,t) as a density time evolution operator; cf. Eq. I|12|) [12J. Poisson's equation is now obtained by eliminating (f> 
from Eqs. lfT3lfTi|l 

a~ — I a -1 — ; I = — 47rZ t ;e(n — fi ti^q) . (16) 



di\ di 

Note that a factor S 2 should appear in the left-hand side if one rescaled Eq. (jlHfl as described above; cf. the last of 
Eqs. jnj. This will be retained for later reference, with respect to the treatment suggested in Ref. |j| (see discussion 
below) . 

In principle, our aim is to solve the system of Eqs. (|12f> to (|15|l or, by eliminating <fi, Eqs. |Q , 1|13|) and 116|) 
for a given initial condition n(£,T = 0) = rc-o(£); an d then make use of the definition in order to invert back to 
the Eulerian arguments of the state moment variables (i.e. density, velocity etc.). However, this abstract scheme is 
definitely not a trivial task to accomplish. 
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III. NONLINEAR DUST-ION ACOUSTIC OSCILLATIONS 

Multiplying Eq. 114fl by «(£,t) and then adding to Eq. I|15|) . one obtains 



dr 



Combining with Eq. (T3J, one obtains 



AirZ^n^o n t) . (17) 



d 2 u 2 „ 

q^2 = - U p,i nU i ( 18 ) 



where u) v .i is the ion plasma frequency (defined above). Despite its apparent simplicity, Eq. I|18|l is neither an ordinary 
differential equation (ODE) - since all variables depend on both time r and space £ - nor a closed evolution equation 
for the mean velocity r): note that the (normalized) background particle density n depends on the potential <fi 
and on the plasma parameters; sec its definition 0. The evolution of the potential t), in turn, involves u(£,t) 
(via the quantity a(£,r)) and the ion density n(£, r). 

Eq. (|18|l suggests that the system performs nonlinear oscillations at a frequency ui — uip^n 1 ^ 2 . Near equilibrium, 
the quantity h is approximately equal to unity and one plainly recovers a linear oscillation at the ion plasma frequency 
^>p,i- Quite unfortunately this apparent simplicity, which might in principle enable one to solve for w(£, r) and then 
obtain {£,t} in terms of {x,t} and vice versa (cf. Davidson's treatment for electron plasma oscillations in Ref. 0; 
also compare to Ref. [8|, setting 7 = therein), is absent in the general (off-equilibrium) case where the plasma 
oscillations described by Eq. (|18|) arc intrinsically nonlinear. 

Since Eq. l|18|l is in general not a closed equation for u, unless the background density n is constant (i.e. independent 
of <j), as in Rcfs. 0,01), one can neither apply standard methods involved in the description of nonlinear oscillators on 
Eq. ltTH|) (cf. Ref. 0), nor reduce the description to a study of Eqs. (|TH1[T7|) (cf. Ref. @), but rather has to retain 
all (or rather five) of the evolution equations derived above, since five inter-dependent dynamical state variables (i.e. 
n, u, E, (j) and a) are involved. This procedure will be exposed in the following Section. 



IV. PERTURB ATIVE NONLINEAR LAGRANGEAN TREATMENT 



Let us consider weakly nonlinear oscillations performed by our system close to (but not at) equilibrium. The basis 
of our study will be the reduced system of equations 

d , \ 
— (an) = 0, 

8u 

dE 

— = (n-n)a, 
aE = -— , 

da du , . 

!Tr=W (19) 



which follow from the Lagrangean Eqs. (|12fl to by scaling over appropriate quantities, as described in ifl"H 
[l3|. This system describes the evolution of the state vector, say S = (a,n,u,E,(f>) (e 5i 5 ), in the Lagrangean 
coordinates defined above. We will consider small deviations from the equilibrium state S = (1, 1, 0, 0, 0) T , by taking 
S = S(°) + eS]/ -* + e 2 S2^°- ) + where e (<C 1) is a smallness parameter. Accordingly, we shall Taylor develop the 
quantity n(4>) near (j> m 0, viz. (f> m e<pi + e 2 (f>2 + in order to express n as 

n PS 1 + Cl cf) + c 2 c/) 2 + c 3 cf) 3 + ... 

= l + ec 1 1 +e 2 (ci(/) 2 + c 2 ^) + e 3 (ci03 + 2c 29 !.i ( /. 2 + C3^) + ..., (20) 

where the coefficients Cj (j = 1,2,...), which are determined from the definition J7J of h, contain all the essential 
dependence on the plasma parameters, e.g. //; making use of e x ~ Y^=o x n /n\, one readily obtains 



d = fi/Zi , c 2 = n/(2Z?) , c 2 = fi/(6Zf) 
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Remember that for /i = 1 (i.e. for vanishing dust) one recovers the expressions for IAWs in e-i plasma. 

Following the standard reductive perturbation technique [l4|, we shall consider the stretched (slow) Lagrangean 
coordinates Z = e(£ — At) , T = e 2 r (where A G Sft will be determined later). The perturbed state of (the j—th 
— j = 1, ...,5 — component of) the state vector is assumed to depend on the fast scales via the carrier phase 

9 = k£ — ujt, while the slow scales enter the argument of the (j—th element's) I— th harmonic amplitude 5^™ , viz. 
Sj^ = X^-oo T) e 1 "^"" 1 ") (where S'j'j'/ = sffl ensures reality). Treating the derivative operators as 

d d d 2 d d d d 

and substituting into the system of evolution equations, one obtains an infinite series in both (perturbation order) e™ 
and (phase harmonic) The standard perturbation procedure now consists in solving in successive orders ~ e™ and 
substituting in subsequent orders. The method involves a tedious calculation which is however straightforward; the 
details of the method are presented e.g. in Ref. 5], so only the essential stepstones need to be provided here. 
The equations obtained for n = I = 1 determine the first harmonics of the perturbation 

n[ 1] =-a[ 1] = (fc 2 /^ 2 )^,"^ =(k/w)ip, e[ 1] = ~ik^ (21) 

where tp denotes the potential correction <j>^ . The cyclic frequency u) obeys the dispersion relation u> 2 = k 2 /(k 2 +sci), 
which exactly recovers, once dimensions are restored, the standard IAW dispersion relation [3j mentioned above. 
Proceeding in the same manner, we obtain the second order quantities, namely the amplitudes of the second 

f2l (2) (2) 

harmonics S 2 and constant ('direct current') terms S , as well as a finite contribution to the first harmonics; 
as expected from similar studies, these three (sets of 5, at each n,l) quantities are found to be proportional to ip 2 , 
|V>| 2 and dip/dZ respectively; the lengthy expressions are omitted here for brevity. The (n = 2, I = 1) equations 
provide the compatibility condition: A = uj(1 — uj 2 )/k = dut/dk; A is therefore the group velocity v g (k) = u'{k) at 
which the wave envelope propagates. It turns out that v g decreases with increasing wave number fc; nevertheless, it 
always remains positive. 

In order ~ e 3 , the equations for / = 1 yield an explicit compatibility condition in the form of a nonlinear Schrodinger- 
type equation (NLSE) 

*||+Pg + Qh/O = 0. (22) 

Recall that ip = <j>± denotes the amplitude of the first-order electric potential perturbation. The 'slow' variables 
{Z,T} were defined above. 

The dispersion coefficient P is related to the curvature of the dispersion curve as P — w"(fc)/2 = — 3w 3 (l — 
w 2 )/(2fc 2 ). One may easily check that P is negative (for all values of k). 

The nonlinearity coefficient Q is due to carrier wave self-interaction. It is given by the expression 



W 12 fc 4 Zf 



3Zfk 6 - 3(a* + 4)Z 2 fc 4 + 3(1 - 2fi - h[i 2 )Zik 2 - ^(3^ - if 



(23) 



where the coefficients 01,2,3 were defined above. 
For low wavenumber k, Q goes to —00 as 



Q 



(3/i-l)V /2 1 



12 Zf 12 



A. Modulational stability analysis 

According to the standard analysis , we can linearize around the plane wave solution of the NLSE (1221 -0 = 
^ e iQM 2 T _|_ c c : [c.c: complex conjugate) - notice the amplitude dependence of the frequency shift Auj = e 2 Q\ip\ 2 
- by setting ?p — -0o + £ $1 j an d then assuming the perturbation ipi to be of the form: ipi — e l ^ k ^ UT ^ + cc 

Substituting into l)22[l. one thus readily obtains Cj 2 = P 2 k 2 (^k 2 — 2(Q / P)\^Ji, \ 2 ^j ■ The wave will thus be stable 
(V fc) if the product PQ is negative. However, for positive PQ > 0, instability sets in for wavenumbers below a critical 
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FIG. 1: A heuristic representation of wave packets modulated by solutions of the NLS equation. These envelope excitations 
are of the: (a, b) bright type (PQ > 0, pulses); (c) dark type, (d) gray type (PQ < 0, voids). Notice that the amplitude never 
reaches zero in (d). 



value k cr — y / 2Q/P |V>i,o|, i- e - f° r wavelengths above a threshold A cr = 2ir/k cr ; defining the instability growth rate 
cr = \Imu}(k)\, we see that it reaches its maximum value for k = k cr /^/2, viz. 

Vrnax = \ Imt2} \k='k ar / V2 = \Q\ 1^1,0 P ■ 

We see that the instability condition depends only on the sign of the product PQ, which may be studied numerically, 
relying on the expressions derived above. 



B. Finite amplitude nonlinear excitations 



The NLSE IL'21) is long known to possess distinct types of localized constant profile (solitary wave) solutions, 
depending on the sign of the product PQ 0, ITEl H(| . Remember that this equation here describes the evolution of the 
wave's envelope, so these solutions represent slowly varying localized envelope structures, confining the (fast) carrier 
wave. The analytic form of these excitation can be found in the literature (see e.g. in [j| for a brief review) and need 
not be derived here in detail. Let us however briefly summarize those results. 

Following Ref. ^(|, we may seek a solution of Eq. (|22[l in the form ?/;(£, t) = p(Z, T) e l0 (£' T ) + c.c, where p, a 
are real variables which are determined by substituting into the NLSE and separating real and imaginary parts. The 
different types of solution thus obtained are summarized in the following. 

For PQ > we find the (bright) envelope soliton 



p = ±po sech 



Z 



L 



(24) 



which represents a localized pulse travelling at the envelope speed u e and oscillating at a frequency 17 (at rest). The 
pulse width L depends on the maximum amplitude square p as L — (2P/Q) 1 /' 2 / p . Since the product PQ is always 
positive for long wavelengths, as we saw above, this type of excitation will be rather privileged in dusty plasmas. The 
bright-type envelope soliton is depicted in Fig. b. 

For PQ < 0, we obtain the dark envelope soliton (hole) ^(| 



P 

e 



±pi 
i 

2P 



1 — sech 



u,,Z 



1/2 



± p\ tanh 



u P T 



V L ' 



2PQ P i T 



(25) 
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which represents a localized region of negative wave density (shock) travelling at a speed u e 
pulse width depends on the maximum amplitude square p\ via L' — {2\P/Q\) 1 / 2 / p\. 
Finally, still for PQ < 0, one also obtains the gray envelope solitary wave 



P = ±P2 



1 — a sech? 



Z — u e T 
U' 



1/2 



see Fig. Again, the 



(26) 



which also represents a localized region of negative wave density. Comparing to the dark soliton (|25[1 . we note that 
the maximum amplitude p2 is now finite (non-zero) everywhere; see Fig. ^1. The the pulse width of this gray-type 
excitation L" = ^J2\P/Q\/ (a pi) now also depends on an independent parameter a which represents the modulation 
depth (0 < a < 1). The lengthy expressions which determine the phase shift and the parameter a, which are 
omitted here for brevity, can be found in Refs. Il6|. For a = 1, one recovers the dark soliton presented above. 

An important qualitative result to be retained is that the envelope soliton width L and maximum amplitude p 
satisfy hp ~ y 1 P/Q (see above), and thus depend on (the ratio of) the coefficients P and Q; for instance, regions 
with higher values of P (or lower values of Q) will support wider (spatially more extended) localized excitations, for 
a given value of the maximum amplitude. Contrary to the KdV soliton picture, the width of these excitations does 
not depend on their velocity. It does, however, depend on the plasma parameters, e.g. here p. 

The localized envelope excitations presented above represent the slowly varying envelope which confines the (fast) 
carrier space and time oscillations, viz. <fi — ^(AT, Z) cos(fc£— lot) for the electric potential (f> (and analogous expressions 
for the density rii etc.; cf. <|21[) h The qualitative characteristics (width, amplitude) of these excitations, may be 
investigated by a numerical study of the ratio P/Q = T](k; p): recall that its sign determines the type (bright or 
dark) of the excitation, while its (absolute) value determines its width for a given amplitude (and vice versa). In 
Fig. |21 we have depicted the behaviour of 77 as a function of the wavenumber k and the parameter p: higher values 
of p correspond to lower curves. Remember that, for any given wavenumber k, the dust concentration (expressed 
via the value of p) determines the soliton width L (for a given amplitude p; see discussion above) since L ~ r/ 1 ^ 2 / p. 
Therefore, we see that the addition of negative dust generally (p < 1) results to higher values of rj (i.e. wider or higher 
solitons), while positive dust {p > 1) has the opposite effect: it reduces the value of n (leading to narrower or shorter 
solitons). In a rather general manner, bright type solitons (pulses) seem to be rather privileged, since the ratio rj (or 
the product PQ) of the coefficients P and Q is positive in most of the k, p plane of values. One exception seems to be 
very the region of very low values of p (typically below 0.2), which develops a negative tail of r\ for small k (< 0.3 ko)'- 
thus, a very high (> 80 per cent) electron depletion results in pulse destabilization in favour of dark-type excitations 
(Fig. d). Strictly speaking, r\ also becomes negative for very high wave number values (> 2.5 kr>)\ nevertheless, we 
neglect - for rigor - this region from the analysis, in this (long wavelength A) fluid picture (for a weak dust presence, 
short A DIAWs may be quite strongly damped; however, this result may still be interesting for a strong presence of 
dust, when Landau damping is not a significant issue Q). 




1.2 



0.4 0.6 0.8 1 1.2 1.4 
wave number k/k D 



FIG. 2: The ratio n = P/Q of the coefficients in the NLSE l|22|l is depicted versus the wave number k (normalized over fcn), 
for several values of the dust parameter p; in descending order (from top to bottom): 0.8, 0.9, 1.0, 1.1, 1.2. 
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RELATION TO PREVIOUS WORKS: AN APPROXIMATE NONLINEAR LAGRANGEAN 

TREATMENT 



By combining the Lagrangean system of Eqs. Ijl2|l to (|16J) . one obtains the (reduced) evolution equation 

d 2 ±d_ 
dr 2 \n J no <9£ 

where the function g(x) [defined in Eq. JSJ] is evaluated at 

d 2 I I 



n d(, 



(27) 



w(n) = n 



1 



dr 2 



Note that the ion density n has been scaled by its equilibrium value n^.o, to be distinguished from the initial condition 
n o = n(£,r = 0). 

Despite its complex form, the nonlinear evolution equation l|27[) can be solved exactly by considering different 
special cases, as regards the order of magnitude of the dispersion-related parameter i5. This treatment, based on Ref. 
, will only be briefly summarized here, for the sake of reference. 

First, one may consider very short scale variations, i.e. L <C Xd (or S 1). This amounts to neglecting collective 
effects, so oscillatory motion within a Debye sphere is essentially decoupled from neighboring ones. By considering 



w{n) w — S 2 nd 2 (l/nj/dr 2 and cj) ~ (i.e. h ~ 1), one may combine Eqs. I|15|l and (|27[) into 



dr 2 



\7, 1 



= 0. 



(28) 



which, imposing the initial condition n(£,0) = ?"io(£), yields the solution 

no(0 



(29) 



Note that if the system is initially at equilibrium, viz. no(£) = ni,0i then it remains so at all times r > 0. Now, 
may go back to Eq. (|12|l and solve for a(£, r), which in turn immediately provides the mean fluid velocity u 



u(4,t) = Up iSmuip iT / 1 d4 

ri it0 J 



via (|lt)|l . and then E(£,t) and 0(£,r). Finally, the variable transformation 10 may now be inverted, immediately 
providing the Eulerian position x in terms of £ and r. We shall not go into further details regarding this procedure, 
which is essentially analogue (yet not identical) to Davidson's treatment of electron plasma oscillations. 

Quite interestingly, upon neglecting the dispersive effects, i.e. setting S — 0, Eq. Ij27(l may be solved by separation 
of variables, and thus shown to possess a nonlinear special solution in the form of a product, say n(£, r) = ni(^)n2(r) 
fl3| . This calculation was put forward in Ref. (where the study of IAW - in a single electron temperature plasma - 
was argued to rely on an equation quasi- identical to Eq. (|27|l '). However, the solution thus obtained relies on doubtful 
physical grounds, since the assumption 5 ~ 0, which amounts to remaining close to equilibrium - cf. the last of Eqs. 

, implies an infinite space scale L (recall the definition of 5) , contrary to the very nature of the (localized) nonlinear 
excitation itself. Rather not surprisingly, this solution was shown in Ref. to decay fast in time, in both Eulerian 
and Lagrangean coordinates. Therefore, we shall not pursue this analysis any further. 



VI. DISCUSSION AND CONCLUSIONS 



We have studied the nonlinear propagation of dust ion acoustic waves propagating in a dusty plasma. By employing 
a Lagrangean formalism, we have investigated the modulational stability of the amplitude of the propagating dust 
ion acoustic oscillations and have shown that these electrostatic waves may become unstable, due to self interaction 
of the carrier wave. This instability may either lead to wave collapse or to wave energy localization, in the form of 
propagating localized envelope structures. We have provided an exact set of analytical expressions for these localized 
excitations. 

This study complements similar investigations which relied on an Eulerian formulation of the dusty plasma fluid 
model 0. In fact, the Lagrangean picture provides a strongly modified nonlinear stability profile for the wave ampli- 
tude, with respect to the previous (Eulerian) description; this was intuitively expected, since the passing to Lagrangean 
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variables involves an inherently nonlinear transformation, which inevitably modifies the nonlinear evolution profile of 
the system described. However, the general qualitative result remains in tact: the dust ion acoustic-type electrostatic 
plasma waves may propagate in the form of localized envelope excitations, which are formed as a result of the mutual 
balance between dispersion and nonlinearity in the plasma fluid. More sophisticated descriptions, incorporating e.g. 
thermal or collisional effects, may be elaborated in order to refine the parameter range of the problem, and may be 
reported later. 
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